Evolution of turbulent spots in a parallel shear flow 
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The evolution of turbulent spots in a parallel shear flow is studied by means of full three- 
dimensional numerical simulations. The flow is bounded by free surfaces and driven by a volume 
force. Three regions in the spanwise spot cross-section can be identified: a turbulent interior, an 
interface layer with prominent streamwise streaks and vortices and a laminar exterior region with 
a large scale flow induced by the presence of the spot. The lift-up of streamwise streaks which is 
caused by non-normal amplification is clearly detected in the region adjacent to the spot interface. 
The spot can be characterized by an exponentially decaying front that moves with a speed different 
from that of the cross-stream outflow or the spanwise phase velocity of the streamwise roll pattern. 
Growth of the spots seems to be intimately connected to the large scale outside flow, for a turbulent 
ribbon extending across the box in downstream direction does not show the large scale flow and 
does not grow. Quantitatively, the large scale flow induces a linear instability in the neighborhood 
of the spot, but the associated front velocity is too small to explain the spot spreading. 
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I. INTRODUCTION 

The transition to turbulence in spatially extended sys- 
tems does not necessarily take place in all points simulta- 
neously but can be preceeded by the formation of local- 
ized structures that grow to eventually fill space. Already 
the first experiments by Reynolds in pipe flow showed the 
formation of turbulent spots and slugs 0J^] . In Couette- 
Taylor flow between counterrotating cylinders turbulence 
can be confined to propagating spirals Localized 
turbulence has also been observed in plane Couette flow 
H where the fraction of space filled with turbulent flow 
has been used as a measure to define the transition to 
turbulence ||-||]. Besides these transitional phenomena 
localized turbulent spots can also be observed in high 
Reynolds number boundary layers j9|. 

It is tempting to connect both the localization of the 
spots, i.e. the coexistence of a laminar and turbulent 
phase of the shear flow, and the propagation of the 
sharp boundaries, i.e. a front-like structure, to phenom- 
ena studied in considerable detail within amplitude mod- 
els [0. Indeed, some models show qualitatively similar 
behavior. There are, however, several problems that raise 
questions about the applicability of such models. For in- 
stance, they are not derived from the Navier-Stokes equa- 
tion and the extend to which they reflect the hydrody- 
namical processes and interactions remains open. Fur- 
thermore, amplitude equations work best if they can be 
applied in a situation of linear instability and small am- 
plitudes jll| , such as the onset of Rayleigh-Benard con- 
vection near the critical point |fl2|| . But many of the tur- 
bulent spots arise in shear flows that are linearly stable, 
at least in the Reynolds number region of interest here. 
Such behavior can be captured in higher-order Ginzburg- 
Landau models [fl3|| , but the required large amplitudes 
complicate a quantitative comparison. Moreover, inves- 
tigations of plane Couette flow show that the turbulent 
state is not stable but can decay spontaneously for lower 
Reynolds number values Q-Q. 



It is our aim here to analyze the evolution of turbulent 
spots in parallel shear flows, in particular their spanwise 
spreading. Our flow has free-slip boundary conditions 
and is driven by a volume force. Despite the change in 
boundary conditions we observe features similar to those 
in experiments on plane Couette flow with rigid bound- 
ary conditions and a linear shear profile: this supports 
the expectation that there are perhaps universal aspects. 
The model is moreover well suited for high resolution di- 
rect numerical simulations with a Fourier-pseudospectral 
method and allows for a detailed investigation of the dy- 
namics in the transitional region. In particular, we focus 
on the characterization of the front which separates the 
laminar and the turbulent region, on the mechanism by 
which it propagates, on the Reynolds number dependence 
of the front speed and on the large scale flow in the lam- 
inar surrounding of the spot. As we will discuss in more 
detail in the appropriate sections these aspects comple- 
ment previous numerical and experimental investigations 
in wall bounded shear flows |l7|-p4[ . 

The paper is arranged as follows. After introducing the 
physical model and the numerical procedures in Sec. II 
we discuss in Sec. Ill the hydrodynamics of the spread- 
ing mechanism in some detail. The properties of the tail 
of the envelope, such as spatial decay and spreading ve- 
locity, are discussed and three different regimes of the 
spreading process are identified. In Sec. IV we discuss 
the results and give a brief outlook. 



II. THE MODEL 

The system we consider here is a shear flow between 
parallel free-slip surfaces and driven by a volume force. In 
the streamwise and spanwise direction periodic boundary 
conditions are applied; in the normal direction the nor- 
mal velocity component vanishes in the two bounding 
surfaces. With lengths measured in units of d/2 (half the 
gap width) the periodicities in streamwise and spanwise 
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directions are both 80. The volume force with a sinusoidal 
dependence in normal direction gives rise to a laminar 
profile with velocities ±Uq at the surfaces. The Reynolds 
number is defined as Re — Uod/(2u). In these units the 
incompressible Navier-Stokes equation for a velocity field 
u(x,t) becomes 



- + (u.V)u = -V,+ -V2u- 
V • u = 



(1) 
(2) 



p(x, t) denotes the pressure and f(x, t) the external vol- 
ume force specified below. 




Turbulent Spot 



FIG. 1. Geometry of the flow. The x axis points in the 
streamwise, y in the wall-normal, and z in the spanwise 
direction. The central plane at y = 1 in which the spread- 
ing is analysed is shaded gray. 

Figure 1 shows the Cartesian coordinate system we use, 
with x pointing in the streamwise, y in the wall-normal 
and z in the spanwise directions. The fluid volume is con- 
fined to < y < 2 with boundary conditions 



dy 



du 
dy 



at the surfaces and periodic boundary conditions in 
downstream and spanwise directions. 

The shear flow is driven by a volume force f = 
7r 2 /(4i?e) cos(7rj//2) e x acting in the x-direction, which 
in the laminar regime sustains a flow Uo = cos(7rj//2) e 2 ,. 
The velocity field u(x, t) is decomposed into this laminar 
flow Uo and a turbulent part v(x, t). As shown already by 
Tollmicn, Uo is linearly stable, thus demonstrating that 
Fj0rtoft's theorem is a necessary but not sufficient con- 
dition for the transition to turbulence [^5). Nevertheless, 
for sufficiently large driving the flow shows a transition 
to turbulence. 

The free-slip or stress-free boundary conditions have 
the advantage that the flow can be represented com- 
pletely by Fourier modes so that robust pseudospectral 
techniques based on a 2/3-rule dealiazing and an adap- 
tive Rungc-Kutta scheme for advancing in time can be 
used [^6],£7| • To account for boundary conditions the flow 
is represented by the Fourier sums 

u x (x, t) = u xk (t) cos(k y y) exp[i(k x x + k z z)\ , (4) 

k 

%(x, t) = 2J %k(*) sm(k y y) exp[i(k x x + k z z)} , (5) 

k 

u z (x, t) = } j Uzk(t) cos(k y y) exp[i(k x x + k z z)\ . (6) 



with wavenumbers 
k 



7T 

°'2'"' 



NyTT 



k x = o,±^,±2 2 l,... ± ^ 2 JL 

2tt 2tt N z 2tt 

k z = 0,± — ,±2 — ,. ..± — — . 

L, L, 2 L~ 



(7) 
(8) 
(9) 
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In [^8[ and ||, respectively, low-dimensional models for 
the transition to turbulence in plane shear flows with 
stress-free boundary conditions were discussed. Their ba- 
sic flow has the form Uo x (y) ~ sm(7n//2) and is con- 
fined to an interval y € [—1,1]. Both expansions can thus 
be related by a shift in the interval, but we prefer (|J)- 
(^) as it has the more compact representation in sines 
and cosines and is easier to implement numerically. The 
spectral resolution for all runs with evolving spots was 
N x x N y x N z = 256 x 33 x 512. The initial localized 
perturbation is a poloidal vortex of the form 



v(x, f = 0)=Vx Vxylexp[-a|(x - x ) 2 - a 2 (y - y a ) 2 - a 2 z (z - z ) 2 )e 



(10) 



positioned slightly off center in order to avoid spurious 
effects due to accidental symmetries. This initial condi- 
tion is a model for the flow induced in experiments where 
a small transverse jet penetrates the laminar shear profile 
in wall- normal direction [^|-(7|. 
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FIG. 2. Statistical analysis of the transition to turbu- 
lence. Panel (a) shows the number N t of initial con- 
ditions that lead to persistent turbulence for different 
Reynolds numbers and for several spectral resolutions. 
Thirty trajectories were run for every Reynolds num- 
ber. The dashed line is a two parameter least square fit 
N t (Re) = 15 x taah[(Re- Ao)/At] + 15. Panel (b) shows 
the temporal relaxation of the turbulent kinetic energy, 
averaged over all trajectories that became turbulent. The 
double headed arrow marks the interval of temporal av- 
eraging used for Eq. (|f). 

For the analysis in the following sections we need the 
Reynolds number above which a transition to turbulence 
occurs. Because of the free-slip boundary conditions at 
the surface it can be expected to be below the one for 
rigid walls. But as in that case the transition is strongly 
intermittent and the best approach to a definition of a 
critical Reynolds number uses a statistical analysis of run 
time experiments with different initial conditions p9[ . At 



each value of Re we run thirty trajectories starting from 
states with slightly different amplitudes A. The different 
initial conditions were obtained by switching on a field 
( |l0| ) with time-dependent amplitude A(t) = ao sin 2 (7rf/2) 
for t <E [0, 2] , where the factor ao was increased from 1.248 
to 1.326 in steps of 0.0026. Since the aim is to obtain in- 
formation on the bulk properties, we took a smaller box 
with aspect ratio L x : L y : L z = 40 : 2 : 20 (in units of 
half the gap width), and a lower spectral resolution. As 
in the direct numerical simulations of plane Couette flow 
p9| , two types of dynamics can be identified: in one case 
the flow builts up and the energy relaxes with oscillations 
to the turbulent state, i.e. E^ n ^ 0, whereas in the other 
case the state decays towards the laminar profile, so that 
finally Eki n becomes negligible. The number Nt of initial 
conditions that became turbulent vs. Reynolds number 
is shown in panel (a) of Fig. 2. For Re > 200 all samples 
relax to the turbulent state and the turbulent energy in- 
creases linearly with Reynolds number, giving rise to the 
relation (p"3|). As in the case of plane Couette flow, the re- 
laxation to the turbulent state is oscillatory (Fig. 2b). We 
conclude from these studies that more than half the ini- 
tial conditions will become turbulent for Reynolds num- 
bers of 130±5, which, as expected, is lower than the value 
of 320 ± 10 for rigid boundary conditions |29|. Below this 
Reynolds number most turbulence is transient and a tur- 
bulent region can disappear by erosion from within. Thus, 
the kind of spot spreading phenomenon we are interested 
can only occur for Reynolds numbers above this value. 



III. SPREADING OF THE TURBULENT SPOT 

After a short transient of about 5 time units the inital 
perturbation which was localized in diameter to about 4 
in half gap width units develops streamwise streaks and 
vortices and starts to expand. To highlight the turbulent 
deviations about the laminar flow a contour plot of the 
downstream velocity averaged over half a box height, 



V x (x,Z,t) 



v x (x,y,z,t)dy , 



(11) 



is shown in Fig. 3 (left panel). The elongated stream- 
wise streaks with alternating flow direction stand out 
above the background flow. The cut at x = L x /2 = 40 
for the streamwise turbulent velocity itself underlines 
the existence of the streamwise streaks (right panel). In 
the streamwise direction the spot advances more or less 
stochastically, with the unpredictable appearance of tur- 
bulent bursts which are then advected by the laminar 
profile. This causes a strongly fragmented spot interface. 
In the spanwise direction it advances more steadily with 
a regular interface and it is this direction we focus on in 
the following analysis. 
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FIG. 3. Streamwise streaks in an expanding spot for Re = 200 at t = 39. Left: Contours of v x (x, z). The rectangular 
box in the spot center marks the lateral extension of the volume used for the analysis of internal turbulent fluctu- 
ations. Right: Magnification of the wall-normal-spanwise plane at x = L x /2 (marked by arrows in the left panel). 
Corresponding contours of v x (x — L x /2, y, z) are plotted using the same linestyle: gray lines denote negative values, 
black lines positive ones, and the heavy solid line is v x = 0. 



Cross sections of the local turbulent energy v 2 taken 
in the middle of the cell at x — L x /2 and y = 1 for fixed 
time indicate three different flow regimes: a turbulent in- 
terior, a laminar exterior and a narrow transitional region 
with large velocity amplitudes and a rather regular spa- 
tial structure (see Fig. 4). We will discuss these regions 
in turn. 



Turbulent interior and wave propagation at the 
spot interface 



The turbulent fluctuations in the interior were inves- 
tigated in detail for five values of the Reynolds num- 
ber, as shown in Tab. [j]. The fluctuations, in units of 
were obtained by averaging over a box V c of size 



0/ 



IT 2 



lx ■ ly ■ lz = 12 : 2 : 6 (see the left panel of Fig. 3) in the 
center of the spot according to 



k 2 )(*o) 



i 



vf(x, y, z, to) d V for i = x,y,z. 



(12) 
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FIG. 4. Envelopes of the turbulent kinetic energy along 
the spanwise axis for three different Reynolds numbers 
at t — 59 after inducing the perturbation. 
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FIG. 5. Snapshots of the streamwise vorticity uj x (upper row) and the streamwise velocity v x (lower row). The data 
are for Re = 200 at x — L x /2 — 40. Dotted lines denote positive contours and solid lines negative ones. Horizontal 
axis denotes wall-normal direction in all panels. 



For a Reynolds number of about 150, where the spot 
decays, the fluctuations decrease with time. For higher 
Reynolds numbers they first increase until a time of 
about 30 and then stay constant, within statistical fluc- 
tuations. The values quoted in the table are the temporal 
averages taken between t = 30 and 50. As expected, the 
streamwise fluctuations are largest but smaller than the 
fluctuations in the interface region (see Fig. 4). 

The spatial modulations near the interface are due 
to elongated streamwise structures, so-called streamwise 
streaks, which are not stable but travel slowly in span- 
wise direction with a phase velocity v w . The occurence of 
such (oblique) waves was noted in several previous inves- 



tigations, mainly for plane Poiseuille flow, and has been 
connected to an inflectional instability of the combined 
shear and cross flow velocity field . A linear stabil- 

ity analysis of the combined profile for plane Poiseuille 
and plane Couette flow gives a range of critical wavenum- 
bers (see Jl8|,|2Cj and below). 

The spreading of the spot in the spanwise direction is 
further documented in Fig. 5 for Re = 200. The stream- 
wise vorticity component u x = d v v z — d z v y is shown in 
the upper panels. The growth of a pair of new counter- 
rotating vortices (panels No. 2,3, and 4 of Fig. 5) can 
be followed for all eight snapshots. The corresponding 
streamwise streak which is lifted up by the non-normal 
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amplification can also be identified in the cross-section 
of the streamwise velocity v x (see the lower panels) . Al- 
though there is a certain discreteness in the growth of 



new streamwise rolls and pairs of streamwise vortices, 
the front advances steadily in time with out any disrup- 
tions as will be shown later in Sec. [II C. 



Re Vw vf \J (v 2 ) y 7 {vj) y/_ {vf) \J 

150 0.02 0.09 0.19 0.18 0.02 0.05 

200 0.08 0.42 0.41 0.38 0.09 0.12 

250 0.14 0.56 0.50 0.45 0.13 0.19 

300 0.17 0.63 0.61 0.56 0.15 0.19 

350 0.16 0.64 0.57 0.53 0.13 0.19 



TABLE I. Characteristic velocities of the spreading spot in units of Uo for the several Reynolds numbers: the phase velocity 
v w , the front velocity vf, and the root mean square velocities of the three components Vi with i — x,y and z. The phase velocity 
was determined separately for the maxima and minima of v x closest to the boundary and their arithmetic mean is listed. The 
turbulent fluctuations taken over a small volume in the spot center were in addition averaged over time between t — 30 and 50. 
Third column is the root mean square velocity, {{v 2 ) + {v'y) + (u^)) 1 ' 2 . 



In Table | we have included the results for the phase 
velocity v w of the streaks. The wavelength was also de- 
termined as the distance between two streaks with the 
same sign, i.e., the spatial distance between two maxima 
or minima of v x , measured in units of d/2. We found val- 
ues between 2.4 and 3.4 with a tendency towards smaller 
values for higher Re. Velocities were determined by mon- 
itoring the motion of these maxima and minima in v x . 
Only extrema closest to the interface were included and 
measurements were limited to times between t = 32 and 
t = 60 for all five data sets. We find that the phase speed 
increases with Reynolds number. 



B. The large scale flow outside the spot 



Near the spanwise centerline across the spot we find a 
strong outward pointing flow. It varies with height but 
is not compensated by an inflow on any level. Incom- 
pressibility thus demands a compensating inflow in other 
parts of the spot. The magnitude of this large scale flow 
decreases rather rapidly with distance from the spot and 
accounts almost completely for the deviations from the 
laminar profile. In order to highlight the flow pattern we 
show in Fig. 6 the directional field, i.e. v/|v|, where the 
overbar indicates an average in the normal direction. The 
flow has quadrupolar characteristics, with outflow in the 
spanwise direction and inflow in the streamwise direction. 

The outward velocity does not coincide with the phase 
speed of the waves. But small as the large scale flow may 
be, it has profound consequences for the spreading of the 
spot. Consider, for instance, the case of a ribbon spanning 
the periodicity box in streamwise direction but localized 
in spanwise direction: no such quadrupolar flow can form. 
And indeed, the ribbon does not spread! (see the dashed 
line in panel (a) of Fig. 9). 
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FIG. 6. Large scale flow outside the spot. To emphasize 
the topology of the flow field only the direction v/|v|, av- 
eraged in wall- normal direction x, is shown. The domain 
is the lower half of the full integration domain, with the 
center of the spot in the middle of the upper boundary. 
The flow is for a Reynolds number of Re = 200 at a time 
t = 39. 



C. The propagating front 

Outside the spot and into the laminar regime one notes 
a gradual decrease in local turbulent energy v 2 (x, t). 

Quantitatively, the energy density decays exponen- 
tially, as demonstrated in Fig. 7, where segments of the 
turbulent spot envelopes of the lower front (moving to- 
ward z = 0) between t = 34 and t = 50 are plotted. The 
rate of spatial decay A grows slowly for Reynolds num- 
bers Re > 200. Fitting to each of 11 snapshots separately 
an exponential profile ~ exp(— z/X) gives decay rates A 
between 0.9 and 1.2. In Fig. 8 the arithmetic mean and 
the corresponding error bars are shown for Re > 175. 
For smaller Reynolds numbers deviations from an expo- 
nential envelope are larger, resulting in larger variations 
and uncertainties in the spatial decay rates. This expo- 
nentially decaying envelope is not much influenced by 
the turbulent fluctuations inside the spot and advances 
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more or less steadily into the laminar region. Monitor- 
ing the position of a certain turbulent intensity thresh- 
old allows to extract a velocity vp that is independent 
of the selected threshold along the tail. However, since 
the turbulent intensity increases with Reynolds number 
it is advisable for numerical reasons to also adjust the 
threshold. We choose to increase the level linearly with 
Re, 

v 2 = 7.5 ■ 1CT 3 x [01 + a 2 (Re - Re )] , (13) 



so 



£ 40 




Re=200 



where the coefficients a\ = 0.143, a 2 — 2.6 ■ 10" 4 , and 
Reo — 200 follow from a linear fit to the mean kinetic 
energy Ek in {t) = (1/2V) / y v(x, t) 2 dV as a function of 
Re. This variation of mean kinetic energy was determined 
alongside with the statistical analysis needed to deter- 
mine the Reynolds number for the transition to the tur- 
bulent state (see section ||) . 
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FIG. 7. The exponential decay towards the laminar pro- 
file outside the turbulent spot at Re = 200. Shown are 
the profiles of the turbulent kinetic energy along the lower 
half of the center line on a semi-logarithmic scale. Curves 
for different times in [34, 50] are vertically displaced and 
overlayed. The lines across the plot indicate the motion 
of the levels at which the front velocity was determined. 
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FIG. 8. Exponential decay rate A of the envelope. The 
data are mean values taken from front envelopes as shown 
in Fig. 7. 
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FIG. 9. Reynolds number dependence of the front veloc- 
ity, (a): Position of the extracted level v 2 of the envelope 
along the spanwise centerline (x = 40, y = 0, z) vs. time 
t for both spanwise fronts. The outer most curves for 
the lower and upper front are taken at Re — 400, the 
inner ones at Re = 125. The dashed lines correspond 
to a run at Re = 200 with a turbulent ribbon that ex- 
tended across the box in downstream direction, (b): Cor- 
responding front velocities vf (absolute values) extracted 
within the shaded time interval in the panel (a). The inset 
highlights the crossover to the shrinking regime around 
Re ~ 135. 
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The front speed thus determined is shown in panel (b) 
of Fig. 9. It increases for Reynolds numbers between 135 
and 200 and saturates for Reynolds numbers above 300. 
Both regimes were also found in plane Couette flow ex- 
periments with rigid walls. Dauchot and Daviaud ]22] ] 
observed an increasing spreading rate for Re between 
370 and 450. Tillmark and Alfredsson found a constant 
spreading rate for Re > 500 [|| . In our simulations we can 
cover both ranges. Note that for Re < 135 the velocity is 
negative, i.e. the spot shrinks, rather than expands. This 
regime is difficult to detect numerically since at these 
Reynolds numbers the turbulence is not very stable and 
can decay spontaneously. This corresponds to an erosion 
of the spot from the inside. However, inspection of the 
flow field shows that in the time interval followed here 
the velocity given is connected with a retreating front 
and not an eroding spot (see also the innermost curves 
for the smallest Reynol ds nu mbers in Fig. 9a). As al- 
ready mentioned in Sec. Ill A the spreading is monotonic 
in time as demonstrated by the curves in Fig. 9a. 

The Reynolds number at which the spots start to grow 
was used by Lundbladh and Johansson |19f | as a definition 
for the threshold to turbulent behavior. Here we find that 
the spot starts to grow for Re ~ 135, a value compati- 
ble with the one determined by the statistical analysis in 
section |l[ so that the results are consistent with Jl9| . 



D. The propagation mechanism 

By now we have identified three velocities near the 
boundary of the turbulent spot: the front velocity Vf 
with which the exponentially decaying part of the enve- 
lope propagates, the velocity of the spanwise outflow U z 
and the phase velocity v w of the streamwise roll pattern. 



The value for U z depends on the distance from the spot, 
but even when calculated at the position of the first max- 
imum of the turbulent fluctuations (see e.g. Fig. 4), the 
value is smaller than the front speed vp. Figure 5 shows 
that the outwards travelling roll patterns are eventually 
overtaken by the turbulent interior since vf is larger than 
V w . 

In the absence of a linear instability of the laminar 
shear flow two possibilities for the growth mechanism 
have to be considered: a linear instability induced by 
the cross flow P,|T7|,|l8|,p2| and non-normal amplification 



The combination of basic profile and cross flow defines 
a rather steady laminar flow with stability characteristics 
different from those of the laminar profile. In the case of 
plane Poiseuille flow, Henningson could show that the 
combined flow is linearly unstable, but the front veloci- 
ties deduced from this instability were smaller than the 
observed spreading velocity Jljjjl^] ■ A complete stability 
analysis would have to include the full profile of the cross 
flow. An estimate of the expected spot growth rates may 
be based on a local approximation, where the values of 
the cross flow are kept fixed. We thus determine for each 
point zq along the spanwise centerline the cross-flow 

U z {y,z a ) = U z (0 1 z a ) + - [1 - cos(Try)] , 

(14) 

where U z (0,z o ) ~ U z (2,z ), and analyze the stability 
against perturbations 

v y (x, y, z, t) = v y (y) exp[i(k x x + k z z - wt)] , (15) 

where ix> = uj r + ie is a complex frequency. This leads to 
the Orr-Sommerfeld equation 



iu(D 2 



k 2 )v y + {ik x U ' x + ik z U z )v v + Re-\D 2 - k 2 ) 2 v y - {ik x U 0x + ik z U z ){D 2 - k 2 )v y = , 



(16) 



with boundary conditions 

v y (v) = % (y) = at y = and 



(17) 



Here, k 2 = k 2 + k 2 . The primes on the basic profiles 
U and D denote derivatives with respect to the wall- 
normal coordinate y. The local approximation is now 
reflected in the fact that the perturbations can have a 
z-dependence, but the basic profile around which the 
perturbations are analyzed does not. With the form of 
the cross flow profile and a Fourier expansion for v y (y) 
the Orr-Sommerfeld equation can be solved algebraically. 
The maximal growth rates e = max(e) thus obtained are 
shown in Fig. 10 for points along the spanwise half-axis 
and for two values of the Reynolds number. For the lower 
value no linear instability is detected. For Reynolds num- 
bers Re > 200 the local growth rate becomes positive, 
indicating a linear instability. 



Investigations of other front- propagation problems, 
usually within a Ginzburg-Landau model, show that it 
is not only the local instability that determines the front 
speed but that the local curvature in wave number space 
has to be included as well. In the absence of a deriva- 
tion of an amplitude equation in a turbulent medium we 
phenomenologically take the amplitude in the Ginzburg- 
Landau equation to model the envelope of the turbulent 
intensity v 2 (z,t), calculated in the middle of the cell at 
x = L x /2 and y = 1. The front then connects a laminar 
state (A ~ 0) with a turbulent one (A ^ 0) Ql3|. For 
all practical purposes the turbulent state is stable and 
the laminar one, composed of the basic profile and the 
cross flow, shows a linear instability. Thus the simplest 
Ginzburg-Landau model with cubic nonlinearity should 
be appropriate, 



d t A = eA + Dd 2 A~b 3 A 3 



(18) 
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The assumption is that the three parameters e, D > 
and 63 > are real. The marginal stability hypothesis 
then predicts a value for the asymptotic front velocity of 



v* = 2VeD. 



(19) 



When taking the maximum growth rate e of our data and 
evaluating the diffusion coefficient D by a saddle point 
approximation around the maximum of the dispersion 
relation lu = u)(k x ,k z ), we end up with a front velocity 
which is about an order of magnitude smaller than the 
observed one. Corresponding to Eq. ([i"9|), one gets, e.g., 
v* ~ 0.06 for e ~ 0.017, D ~ 0.056, and Re = 300 (see 
Fig. 10). However, at this Reynolds number the front 
moves with vp ~ 0.63. 
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that could arise in such models, is rather simplified and 
can only give an indication of the expected front veloc- 
ity. Nevertheless, it seems to us that even when these 
improvements are included the model cannot account for 
the observed spot spreading rates: first of all, it cannot 
say anything about the dynamics below a Reynolds num- 
ber of about 200, where the spot spreads but where there 
is no instability. And secondly, even above this Reynolds 
number, where there is a linear instability, the calculated 
front speed differs significantly from the observed one. 

Besides the linear instability mode for spot growth 
there is another possibilty, based on the non-normal am- 
plification of perturbations near the spot interface. Be- 
cause of the action of the perturbation on the basic pro- 
file, streamwise vortices need not decay monotonically 
but can first grow on a time scale of order Re to an ampli- 
tude about a factor Re larger than the initial amplitude. 
This so-called lift-up effect |3(| is most likely responsi- 
ble for the occurence of streamwise streaks in turbulent 
shear flows where longitudinal modulations can give rise 
to secondary instabilities and a perpetual non-periodic 
time evolution jl6|,[||]. Indeed we do observe this cyclic 
reproduction of the coherent structures in the interface 
region of the spot as demonstrated in the sequence of 
Fig. 5. The model then is that the turbulent interior of 
the spot induces a small perturbation near the interface 
which will then be amplified and grow turbulent. The sta- 
tistical analysis of section [n] shows that not all perturba- 
tions grow turbulent. Thus, if most of the perturbations 
grow turbulent, the spot will spread, but if most of the 
perturbations decay, it will shrink. The observed coinci- 
dence between the critical Reynolds number for transi- 
tion and the one for spot growth can naturally be ex- 
plained in this picture. However, we do not see how to 
derive other quantitative conclusions from this model. In 
particular, the front speed would be given as the quotient 
of the width of the rolls generated (this can be read off 
rather accurately from the frames) and the time a per- 
turbation needs to grow turbulent. The latter depends 
on the amplitude of the initial seed, the threshold for 
the transition to turbulence and the amplification rate, 
neither of which seems accessible to independent deter- 
mination. 



FIG. 10. Maximum growth rates e of plane waves per- 
turbing the mean flow combined with the cross-stream 
outflow along the spanwise axis. The sets for two differ- 
ent Reynolds numbers are compared at t = 39. Thick 
solid lines denote log e (w 2 ) and thick dotted lines denote 
maximum growth rate e. Additionally, for Re — 300 the 
diffusion coefficient D, which follows from the dispersion 
relation, is plotted for the gray shaded range of zq values. 

This analysis of the front propagation mechanism, with 
a local approximation, only cubic terms in the amplitude 
model and without a discussion of slow transients |341 



IV. FINAL REMARKS 



Our analysis of a shear flow with free slip boundary 
conditions on two parallel surfaces has revealed many 
similarities to plane Couette flow between rigid walls. 
As in that case three velocities connected with the spot 
can be identified, the velocities of the advancing front, 
of the outward flow component and of the phase speed 
of (oblique) waves. They differ in value and in Re- 
dependence, and relations between them are unknown. 
The dependence of the front velocity with respect to 
Reynolds number is consistent with experimental find- 
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ings |5|,|22| . Many of the results reported here parallel the 
ones for plane Poiseuille flow. We also find waves and in- 
stabilities in the neighborhood of the spot, but they do 
not lead to quantitative predictions for the front velocity. 

The conclusions we draw from this investigation high- 
light a dilemma: On the one hand side the large scale 
flow outside the spot does not seem to be important: the 
spot grows independent of whether there is a linear in- 
stability of basic flow plus large scale exterior flow or not, 
and if there is a linear instability the derived front speed 
is slower than the observed one. Moreover, the behavior 
expected from the non-normal amplification mechanism 
can explain some aspects of the dynamics, especially for 
lower Reynolds numbers. On the other hand, if the out- 
flow is suppressed, as in the case of a turbulent ribbon 
that spans the cell in streamwise direction, no growth is 
observed despite the random initialization of seeds near 
the spot interface. Perhaps this can explain the observed 
stop of growth in some experiments . So it seems that 
spot growth is a subtle interplay between local features 
(e.g. non-normal amplification) and global features (such 
as the external flow). The connection between both and 
a quantitative estimate of the growth velocity remain a 
major puzzle in the dynamics of turbulent spots in par- 
allel shear flows. 
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